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I. INTRODUCTION 

A challenging problem in nuclear physics is to calculate 
nuclear properties microscopically, using realistic nuclear 
forces. For light systems {A = 2,3), the Schrodinger 
equation can be solved with a high degree of accuracy and 
both ground state properties and reaction cross sections 
have been calculated 0, II H II H II 0] ■ 

For heavier systems, the Green's function Monte Carlo 
(GFMC) method IMIal as well as the no-core shell model 
(NCSM) approach [l^ have been successfully applied for 
the ah initio description of nuclear properties, using re- 
alistic nucleon-nucleon (NN) and three- nucleon forces. 
GFMC calculations, while more accurate, are limited to 
masses up to A = 12, while the NCSM can handle a larger 
range of masses, up to A = 16 and beyond [TllIT^ . 

However, for A> 3 a big obstacle is encountered if one 
wants to calculate reaction observables involving states 
in the continuum, because of the enormous difficulties in 
calculating many-body scattering states. 

Such difficulties can be avoided using integral trans- 
forms, which reduce the continuum problem to a bound- 
state-like problem El El El, so that only bound state 
techniques are required. While not all i nteg ral transforms 
are appropriate for a practical use |la . Il7| (the inversion 
of the transforms is an "ill posed proble m" |18l ). the in- 
tegral transform with a Lorentz kernel jl9l | appears to 
be the practical tool for such calculations. In fact, it 



*On leave from the National Institute for Physics and Nu- 
clear Engineering "Horia Hulubei", Bucharest, Romania, email: 
istetaemail . arizona.edu 
^email : sof iaOphysics . arizona.edu 
■l-email: s.bacca{Bgsi .de 



has allowed the calculation of electromagnetic reaction 
cross sections beyond break-up thresholds of nuclei from 
A — i I23 to A=7 I2II. Such calculations have also been 
possible thanks to the use of the effective interaction hy- 
perspherical harmonic (EIHH) [2^ l23l | approach, a very 
accurate bound state technique, which, similar to the 
NCSM approach, uses the concept of an effective interac- 
tion to speed up the convergence of the basis expansions. 
In particular, by means of the Lorentz integral transform 
(LIT) method combined with the EIHH technique, total 
photoabsorption cross sections of six- [3 IHl a-nd seven- 
body nuclei [2II with semirealistic NN potentials have 
been calculated. In a very recent work [23 the "^He total 
photoabsorption cross section could even be calculated 
with a realistic nuclear force (two- and three-body po- 
tential). Moreover, a similar formalism has been recently 
applied to describe exclusive electromagnetic processes of 
the four-body nuclear system with a semirealistic poten- 
tial mmii. 

While the EIHH and NCSM approaches are rather sim- 
ilar, only the latter has made use of realistic interactions 
in calculations with A > 4 . Indeed, NCSM has the 
advantage that one can use an equivalent Slater deter- 
minant basis, allowing description of a larger range of 
masses. Therefore, it is of great interest to investigate 
the possibility of applying the NCSM to the LIT equa- 
tions. In fact, if one could obtain reliable results in this 
way, the possibility to study also reactions on heavier nu- 
clei by means of realistic ab initio approaches would open 
up. 

The purpose of the present work is to investigate the 
applicability of the NCSM approach to the solutions of 
the bound-state equations required by the LIT method. 
The similarities between the EIHH and the NCSM makes 
the task straightforward, in principle. On the other hand, 
the practical implementation of the method, especially 
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concerning the problems of convergence, might lead to 
difficulties. 

In the following, we present the results of a test con- 
sisting in calculating, within the EIHH and NCSM ap- 
proaches, the LIT of the four-body response functions to 
two different excitation operators. Because of the con- 
vergence properties, the input interaction is the simple 
Minnesota (MN) [S^ potential. Two operators of differ- 
ent nature and range have been chosen, i.e., the isovector 
dipole and isoscalar quadrupole. 

The paper is organized as follows. In Sec. m the theo- 
retical approaches are discussed, that is the LIT as well as 
the NCSM and EIHII methods, respectively. Results are 
presented in Sec. IIIII while conclusions and perspectives 
are discussed in Section Hvl 



in terms of the Lanczos coefficients a„ and 6„, where 



II. THEORETICAL OVERVIEW 

A. The LIT approach 

The inclusive cross sections of reactions induced by 
perturbative external probes are generally written in 
terms of the so-called response functions defined as: 



rf*/ (^/lOI^-o) SiEf-Eo-uj), (1) 



where lo represents the energy transferred by the probe, 
and O the excitation operator. Wave functions and en- 
ergies of the ground and final states of the perturbed 
system are denoted b y | v^ n/f) and Eq/j, respectively. 

In the LIT method [13 one obtains R{lli) after the in- 
version of an integral transform with a Lorentzian kernel 



Riio) 



■ (2) 



The state 4* is the unique solution of the inhomogeneous 
"Schrodinger-like" equation 



{H ^Eo-(7R + i(Ti)\^)^0\^^) 



(3) 



Because of the presence of an imaginary part tr/ in Eq. Q 
and the fact that the right-hand side of this same equa- 
tion is localized, one has an asymptotic boundary con- 
dition similar to a bound state. Thus, one can apply 
bound-state techniques for its solution, and, in particu- 
lar, expansions over basis sets of localized functions. In 
the present paper we solve Eq. ^ using both the NCSM 
and EIHII methods, which will be briefly described be- 
low. 

In both cases we evaluate the LIT by calculating the 
quantity L{a) = (^'1^') directly [sj via the Lanczos al- 
gorithm. Hence, one ffiids that the LIT can be written 
as a continuous fraction 



L[a) = ^ Im 

"I [z - ao) 



1 



(4) 



f — o'R + i<^i and z = Eq + cr . The quantity 
Mo = (*o|OtO|*o) 



(5) 



is the zero-th moment of the distribution Riuj) and, 
therefore, equivalent to the total strength induced by 
the excitation operator (non-energy weighted sum rule). 
A parameter-free smooth response function R{uj) is ob- 
tained by inversion of the integral transform of Eq. Q 
(for inversion methods see Refs. H2 and [13 )■ 

Here we would like to note that similar applications 
have already been used in the past [3 US for evaluating 
response functions by means of shell-model calculations 
in terms of the so-called "Lanczos response" (LR). In 
these cases, LR is 



LR{uj) 



n=l 



SiE„ -Eo-uj) 



(6) 



where Af is the number of Lanczos iterations, while £"„ 
and Sn are the eigenvalues and transition strengths ob- 
tained with the Lanczos algorithm, respectively. Starting 
with Eq. one usually replaces the 6 function by a res- 
olution function (Lorentzian or Gaussian) with a width 
parameter appropriate to the experiment (a typical value 
is ai = 0.25 MeV) js^, or one evaluates the running in- 
tegral 



Eth 



LR(ijj') duj' , 



(7) 



(z-ai)- 



where Eth is the break-up threshold energy of the system. 
While the LR has the same underlying physical content 
as L((Tfl, (T/), since in both approaches one starts from 
the poles and the residues of the Green's function, it is 
evident that the information is processed in a different 
way. In fact, in the LIT method Eq. (0J is regarded as 
an integral transform of the response function, which is 
then recovered via numerical inversion. In this approach, 
one typically chooses larger values of ct/ ( 10 — 20 MeV), 
of the order of the strunctures expected in the response 
function, such that L(aR,ai) can be calculated with the 
sufficient numerical precision that allows a stable inver- 
sion. After the inversion one recovers a smooth response, 
which is independent of aj in a large range of values for 
this parameter. 



B. The NCSM and the EIHH approaches 

The EIHH and NCSM approaches are spectral resolu- 
tion methods, where one performs an expansion of the 
Schrodinger wave function in terms of a complete set of 
basis states: the hyperspherical-harmonics (HH) func- 
tions in the EIHH case, and the HO basis functions in 
the NCSM case. Obviously, due to computational lim- 
itations, the basis has to be truncated. The truncated 
set of states forms the so-called model or P-space, and 
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the excluded basis functions build its complementary Q- 
space, such that the sum P + Q = \ spans the total 
Hilbert space. It is well-known that for a given NN in- 
teraction the HO and the HH expansions converge rather 
slowly. In order to accelerate the convergence pattern, 
both methods make use of an effective interaction intro- 
duced by Navratil et. al. [s^lssf . which follows a unitary 
transformation approach. For the finite P-space, the cor- 
responding effective interaction is obtained by means of 
a unitary transformation [3^|43,0,E3- In general, the 
exact effective interaction for an A-hody system will have 
irreducible A-body matrix elements, even if the bare in- 
teraction is only two-body, and its construction requires 
the exact solution of the initial problem. Since the solu- 
tion for the A-hody problem is the primary goal of the 
approach, an exact calculation of the effective interaction 
is not practical, and one has to resort to approximations. 
Thus, the effective interaction is approximated as a sum 
of a-body effective interaction terms (a < A), derived 
from the solution of the a-body problem. This procedure 
is called the cluster approximation, and in practice, we 
use a = 2 or 3. The effective Hamiltonian is then ob- 
tained by replacing the bare interaction in the A-body 
Hamiltonian by the effective one. However, because of 
the approximation introduced in order to compute the 
effective interaction, the resulting eigenvalues of the ef- 
fective Hamiltonian do not exactly reproduce the values 
in the full space, the signature being a dependence of the 
observables on the parameters of the model space. For a 
given cluster approximation, the eigenenergies converge 
to the bare ones by increasing the size of the model space. 
And in the limit P — > 1 this procedure becomes exact 
|43| . Alternately, one can obtain convergence in a fixed 
model space by increasing the cluster size, as in this case 
the approximate effective interaction converges to the ex- 
act one. 

In the NCSM, one works in an HO basis, and the P- 
space is spanned by states with the total number of HO 
quanta N < A^max- For the purpose of this paper, we 
neglect three-body forces and use local interactions, al- 
though the non-local interactions do not present an im- 
pediment 01- Therefore, the intrinsic Hamiltonian de- 
scribing a system of A nucleons is simply 



H 



A 

E 

i<j=i 



2m 
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Thus, the A-body Hamiltonian takes the form 

A 
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1=1 

A 
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i<j=i 
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2m 
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mfl^r. 



mQ 
~2A 



(10) 



where il is the HO parameter. As the NN potential de- 
pends on the relative coordinates, the added center-of- 
mass HO term has no influence on the internal motion in 
the full space. 

In the present NCSM calculations, we use the three- 
body cluster approximation for the effective interaction. 
In order to derive it, one solves the three-body problem 
for the cluster Hamiltonian 



E 

1=1 

3 

E 

i<j=i 
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2m 



1 



mVl^r,^ 



mVr 



(11) 



and determines the three-body unitary transformation 
from the condition that the transformed Hamilto- 



(12) 



preserves the solutions of Eq. Hll|) in a subspace P*-^-* of 
the full three-body Hilbert space. This approximation 
introduces dependence on the oscillator parameter fl so 
that one has to search for a range of r2-values that mini- 
mizes such a dependence. The effective Hamiltonian for 
the A-body system in the three-body cluster approxima- 
tion is then defined as 



eff 



A 

E 

i=l 



P£_ 

2m 



1 



A 



A 



E V, 

i<j<k=l 



(3) 
ijk ' 



with 



V(3) = 7^(3) _ J2 



Pi 

2m 



1 



2^ 2 



(13) 



(14) 



with Pi and being the momentum and the position 
vectors of the ith particle, respectively; V being the NN 
potential; and m the nucleon mass. In order to facilitate 
the use of the convenient HO basis for evaluating the 
effective interaction, Eq. (jS)) is modified by adding an HO 
center-of-mass Hamiltonian P^/{2Am) + Amn'^R'^/2, 



A A 
P = Y.A and R-jJ2 



(9) 



Note that, because of the decoupling condition 

Q(3)^(3)p(3) ^ ^ 



(15) 



the three-body effective interaction is energy- 
independent. A more detailed description of the 
derivation of the operator X^'^^ can be found in 
Refs. lillillillil. 

In a consistent approach, one should use effective tran- 
sition operators, derived by means of the same transfor- 
mation that determines the effective interaction. While 
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the renormalization at the three-cluster level has not 
been hivestigated for general operators, studies of the 
latter in the two-body cluster approximation have shown 
little effect of the renormalization for long-range observ- 
ables ^3 . Although it is conceivable that a renor- 
malization at the three-body cluster level can show some 
improvement for long-range operators in smaller model 
spaces, this task is extremely demanding and not justi- 
fied since we obtain model space independent results in 
the largest model spaces used in this work. 

The present four-nucleon P-space calculations are per- 
formed in a properly antisymmetrized translationally in- 
variant HO basis, as described in 0|, where the inter- 
ested reader can find details on how to compute an one- 
body translationally invariant operator in such a basis. 

In the EIHH method the calculation is per- 

formed with the HH basis and the spaces P and Q are 
defined by means of the hyperspherical quantum number 
K. Thus, the model space includes all the HH functions 
with K < i^niax- The intrinsic Hamiltonian in the hy- 
perspherical coordinates is given by 
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FIG. 1: (Color online) *He ground-state energy for the MN 
potential as a function of the HO excitations allowed in the 
NCSM model space Nmax and of the maximal value of the 
HH grand-angular momentum quantum number Kmax in the 
EIHH expansion. 



III. RESULTS 
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1 

2m 



i<] = l 
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where p is the hyperradius and Ap contains derivatives 
with respect to p only. The grand-angular momentum 
operator is a function of the variables of particles A 
and {A — 1) and of Ka-2i the grand ang ular momentum 
operator of the residual subsystem |50||. Unlike in the 
case of NCSM calculation, in this work we use a two-body 
cluster approximation for EIHH, which already yields an 
excellent convergence pattern. The starting point for the 
construction of the effective interaction is a two-body-likc 
Hamiltonian, 



(17) 



where rA,A-i is the relative coordinate between particles 
A and A — 1, but the total hyperspherical kinetic energy 
is considered. Due to the presence of the hyperradius p, 
which is a collective coordinate, the resulting effective in- 
teraction produces a sort of "medium correction" . More- 
over, since it depends also on the Ka-2 quantum number 
of the A ~ 2 particle subsystem, it is state dependent. 
The sum of these leads to a faster rate of convergence. 
The HH functions arc constructed starting from the rel- 
ative Jacobi coordinates, thus the center of mass motion 
is removed from the very beginning. The Fcrmionic ba- 
sis states are obtained by coupling spin-isospin states, 
i.e., states that belong to irreducible representations of 
the permutation group Sa, with spatial HH states to 
yield completely antisymmetric basis functions. The ba- 
sis is constructed by using the powerful algorithms of 

m mm mil. 



The aim of the present work is to investigate the reli- 
ability of the NCSM approach for the description of in- 
clusive response functions via the LIT method, and not 
to give a realistic prediction of those response functions. 
Therefore, although applications of the NCSM and the 
EIHH have already been performed and published using 
realistic NN forces as well as theoretical three nucleon 
forces, for the purpose of our comparison it is more con- 
venient to employ the semi-realistic MN |33] potential. 
The advantage of the MN potential is that it consists of 
Gaussian-type potentials and has a rather soft core, lead- 
ing to good convergence rate for both the NCSM and the 
EIHH. 

We present results for two operators different in isospin 
nature and range, giving the leading contribution to low- 
energy electromagnetic reactions. These are the isovector 
dipole and the isoscalar quadrupole operators, respec- 
tively: 




i—1 



i—l 



(18) 
(19) 



The NCSM calculations are also carried out for four 
different HO frequencies {hfl = 12, 19,28 and 40 MeV), 
in order to study the dependence of the resulting LIT 
on this parameter and to provide an estimate for the 
theoretical uncertainties. 

While not the goal of the present calculation, we illus- 
trate the convergence patterns for the *He ground-state 
energy in Fig. ^ This convergence is important because 
(i) the ground-state wave function l^fg) enters the source 
term in the LIT equation ||2Jl, and (ii) the total strength 
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TABLE I: Summary of *He ground-states properties obtained 
using the MN |3(ll | potential model without Coulomb inter- 
action. The evaluation of total strength of the isoscalar 
quadrupole transition, Mq{Q), for the NCSM does not take 
into account the hQ = 40 MeV result, which is not yet con- 
verged. 



1.5 





Ea [McV] 


Mo{D) [fm^] 


MoiQ) [fm^] 


EIHH 


-30.779(1) 


0.7883(1) 


8.54(1) 


NCSM 


-30.80(5) 


0.786(6) 


8.56(14) 



Mq depends entirely on I^'o)- A good convergence is 
reached already at K^ax = 10 for EIHH, whereas one 
obtains slower or faster convergence for NCSM, depend- 
ing on the value of the HO frequency CI. In particular, 
one finds a weak dependence on the model space size for 
hn — 19 and 28 MeV. The five best results, which occur 
for Nrnax — 16 and Kmax = 16, all agree within 0.2%. 
Note that the Coulomb interaction is not included, so 
that the ^He g.s. energy obtained in our calculations 
differs from previously published results for the MN po- 
tential (-29.96(1) McV HI). 

In Table m we compare the ^Hc ground state proper- 
ties obtained using the MN two-body interaction in both 
EIHH and NCSM. One finds an excellent agreement for 
the binding energy, as well as dipole and quadrupole total 
strengths, which will be discussed in more detail below. 



A. Dipole response 

We start the discussion with the results for the isovec- 
tor dipole transition. In Fig. |21 we compare the total 
strength Mq, which enters Eq. (j2J), for O = Z?. Wc evalu- 
ate Mq for the dipole operator by means of an expansion 
over basis states with J^T = 1~1. Alternatively, using 
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FIG. 2: (Color online) *He total dipole strength for the MN 
potential as a function of the HO excitations allowed in the 
NCSM model space Nmax and of the maximal value of the 
HH grand-angular momentum quantum number Kmax in the 
EIHH expansion. 
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FIG. 3: (Color online) The LIT [see Eq. @] for the isovector 
dipole transition as a function of an for ai = 20 MeV: con- 
vergence of the EIHH calculation with respect to the maximal 
value of the HH grand-angular momentum quantum numbers 
Kmax /{Kmax + 1) iu the EIHH expansions. 



the non-energy weighted sum rule, one can evaluate Mq 
directly on the ground state as an expectation value of 
two long-range operators (i.e., the mean square charge 



radius 



and the mean square proton-proton radius 



Mq 



Z{Z-l) 



ppi 



(20) 



As expected, the sum rule is satisfied for each model 
space. The rate of convergence of the total strength 
depends upon the HO frequency, as it does for the g.s. 
energy and other observables, and the least dependence 
on the model space is obtained again for Ml = 19 and 
28 McV. Also, the best EIHH and NCSM results {Nmax 
or Kmax = 16) agree within 0.1% for Ml = 12, 19 and 
28 McV, whereas one finds a larger deviation (0.8%) for 
Ml = 40 MeV, not yet completely convergent. The ori- 
gin of this behavior is related to the fall-off of the wave 
function for the HO potential well. As the HO potential 
well becomes steeper and steeper, the wave functions in a 
fixed model space go faster to zero; as a consequence, one 
finds smaller root-mean square radii for the ^He g.s., and 
the expectation value of a long-range operator, like the 
isovector dipole, is more poorly represented. In princi- 
ple, we could improve the convergence of Mq for Ml = 40 
MeV by calculating the expectation value of D'^ D di- 
rectly on the ground state up to Nmax ~ 20. We would 
like to point out that, despite this, the discrepancy never 
exceeds the 1%. 

The overall convergence behavior of the LIT [see 
Eq. ©J is influenced by the convergence of its two com- 
ponents: (i) the total strength Mq, which was discussed 
earlier; (ii) the residual continued fraction of Lanczos co- 
efficients. While the first is governed by the convergence 
of the ground state (non-energy weighted sum rule) , the 
results for the continued fraction have a further depen- 
dence on the J'^T = 1~1 model space. For this reason. 
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FIG. 4: (Color online) The LIT [see Eq. ©J for the isovector 
dipole transition as a function of gr at aj — 20 MeV: con- 
vergence with respect to the HO excitations allowed in the 
NCSM model space for il = 19 MeV (upper panel), and de- 
pendence on the HO frequency Q of the best {Nmax = 16/17) 
NCSM results (lower panel). The solid line represents the 
EIHH resuh for K^ax = 16/17. 



the convergence pattern of the LIT has to be studied as 
a function of the model space for the ground state as well 
as for the LIT state |\I'). Due to the parity difference of 
the two states, the values of K^ax and Nmax chosen for 
the model spaces are even for the ground state and odd 
for J'^ = 1~1 state, respectively. 

Figure shows that the convergence of the LIT for 
EIHH is fast and accurate over the entire range of aR. 
The same plot for the NCSM results at Hfl = 19 MeV, 
presented in the upper panel of Fig. O reveals a slower 
convergence rate. In particular, we notice that, unlike 
for the EIHH, the high tails of the NCSM results 
show oscillations and a satisfactory quenching of such 
oscillations is obtained only for the biggest model space 
(Nmax = 16/17). Furthermore, one finds a stronger de- 
pendence on the model-space size, even in the peak region 
(see the insets of Figs. O and 0J). 

The influence of the HO frequency on the LIT is shown 
in the lower panel of Fig. 21 In a given model space, for 
smaller HO frequencies one obtains a better sampling of 
the complex-energy continuum, and therefore a more ac- 
curate convergence of the LIT, especially in the low-cr^j 
region; on the other hand, faster convergence is achieved 
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FIG. 5: (Color online) The NCSM {Nmax = 16/17) inclusive 
response to the isovector dipole transition [see Eq. as a 
function of uj for four different values of the HO parameter Q: 
comparison with the EIHH {Kmax ~ 16/17) result. 



when the HO length parameter is chosen close to the size 
of the state to be described, and, therefore, larger h^l val- 
ues arc preferrable for describing the ground state. As 
a consequence, in the particular case of the '*He nucleus, 
frequencies in the range 12MeV < hfl < 28MeV repre- 
sent a good compromise, as shown in Fig^ (lower panel). 
Among the four frequencies adopted to evaluate the LIT 
with the NCSM, HQ = 19 MeV leads to the smoothest 
convergence pattern and to the best agreement with the 
EIHH curve, which appears as a solid line in Fig. 0] 

The different behavior of the two methods with respect 
to the size of the P-space used in the calculation is related 
both to the differences in the shape and asymptotic be- 
havior of the HO and HH base states [58ll59| . and, more 
importantly, to the actual number of states, included in 
the expansion. We would like to point out that, due to 
the extra flexibility of the HH basis in the hyperradial 
part of the expansion, for the same value of Nmax and 
Krnax ^ the number of totally antisymmetric basis states 
in the two approaches is significantly different. For ex- 
ample, in the model space defined by Nmax = 16, one has 
2, 775 four-body states for a NCSM calculation, while for 
Kmax = 16, the total number of four-body states in a 
EIHH calculation is 10, 890. Further details on the con- 
nection between the number of basis states in the two 
approaches can be found in Refs. [23 and j60f . 

Figure shows the EIHH and NCSM results for the 
inclusive response to the isovector dipole excitation ob- 
tained by inverting the LIT's of Figs.|21and^ (Note that 
the same results for the response are obtained by invert- 
ing the LIT for ct/ 10 and 15 MeV.) The HO frequency 
value of Ml ~ 19 MeV leads to the best agreement with 
the EIHH response. Indeed, the discrepancy between the 
two numerical results does not exceed 5% in the energy 
intervals immediately close to the disintegration thresh- 
old as well as in the dipole resonance peak region. While 
the resonant peak is equally well described for Ml = 12 
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FIG. 6: (Color online) *He total quadrupole strength for the 
MN potential as a function of the HO excitations allowed 
in the NCSM model space and of the maximal value of the 
HH grand-angular momentum quantum number Kmax in the 
EIHH expansion. 



and 28 MeV, these HO frequencies yield a larger discrep- 
ancy (at most 10% and 15%, respectively) within 3 MeV 
from threshold. As expected from previous observations, 
poorest agreement is found for Ml = 40 MeV, although 
mainly in the low-energy part of the response. In the 
range 60 MeV< w < 80 MeV all four NCSM responses 
agree within the 7% or better of the EIHH result. 



B. Quadrupole response 

We present now the results for the isoscalar quadrupole 
transition. As expected, the total strength Mq for the 
isoscalar quadrupole, shown in Fig. |S1 presents a slower 
convergence rate than the same quantity for the dipole 
operator. The greatest sensitivity to the model-space 
size is obtained again for Ml — 40 MeV. The calcula- 
tion of the total strength for the quadrupole transition 
employs an expansion over an intermediate set of basis 
states with J^T = 2+0. In particular, for the NCSM val- 
ues in the largest model space {Nmax — 16/16), we omit 
the contribution of the transitions to the J^T = 2+0 
states with N^ax = 18, as a full calculation in this case 
would be out of reach. In principle, one could evaluate 
the loss of strength induced by such an approximation 
(introduced due to computational limitations), using the 
non-energy weighted sum rule and calculating the expec- 
tation value of the operator Q'^Q directly on the ground 
state, as was done for the dipole response. However, the 
form of this operator makes the calculation technically 
difRcult. As for the dipole transition, the EIHH calcula- 
tion shows a fast convergence. The agreement with the 
NCSM results is within the 1% for Ml = 12, 19 and 28 
MeV, while hil ~ 40 MeV, for which the convergence is 
not completely reached, shows a 5% discrepancy. Hence, 
we can assume that the transitions to Nmax = 18 can be 
neglected. 



> 



2.0 
1.5 
1.0 
0.5 
0.0 



1 ' 1 ' 1 


' 1 ' 1 ' 1 






■ — hn = 12 MeV 


- 


- / \ 


— - hn = 19 MeV 


- 


- /• \ 


■ — ■ n\l = Zo MeV 






■ ■ ■ ■ nn = 40 MeV 


- 


7 Nma. = 16/16 \ 


EIHH 


— 






_ K,„a. = 10/12 






ai = 20 MeV 






1.1,1 


.1.1.1 





20 



40 



60 80 

an [MeV] 



100 120 



^5.0 



> 



4.0 



.S 3.0 



3 



1.0 



0.0 





' 1 ' 1 ' 1 


1 




■ - Ml = 12 MeV 






- - Ml = 19 MeV 






■ - ■ fin = 28 MeV 






EIHH 






^ MN 




/, 1 . 1 


. 1 rr-i — 





20 30 40 50 60 70 80 

LO [MeV] 

FIG. 7: (Color online) The dependence on the HO frequency 
Q, of (upper panel) the LIT for ai = 20 MeV, and (lower 
panel) the inclusive response for the isoscalar quadrupole 
transition. Note that Nmax = 16/16 is not a complete model 
space (see text for details). The solid lines represent the EIHH 
results for K,nax = 10/12. 



The comparison of the LIT's in the biggest model space 
{Nmax = 16/16, Kmax = 10/12) for the EIHH and the 
four different NCSM calculations is presented in the up- 
per panel of Fig. [7| All the curves, excluding the re- 
sult for hfl = 40 MeV, still far from convergence, show 
a good agreement, and the corresponding inversions are 
compared in the lower panel of Fig.0 The discrepancies 
in the response functions do not exceed the 10% in the 
range 25 MeV < w < 60 MeV, where the curve shows 
a resonant behavior. More delicate are the regions close 
to threshold and high-energy tail, where the response is 
small and the inversion procedure is more sensitive to 
the numerical error present in the LIT. One finds the 
best agreement for hfl = 12 MeV. 



IV. CONCLUSIONS AND OUTLOOK 

We have presented the results of an ab initio calcula- 
tion of the '''He response functions to the isovector dipole 
and isoscalar quadrupole excitations, obtained by means 
of the LIT method, within both the NCSM and the EIHH 
approaches. As NN interaction, we have used the semire- 
alistic MN potential model. The aim has been to inves- 
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tigate the reliability of the NCSM to the description of 
inclusive response functions via the the integral trans- 
form method with a Lorentzian kernel. The chosen tran- 
sitions have allowed us to perform the test for two op- 
erators different in isospin nature and range. In both 
cases our model study has shown that the NCSM can be 
successfully applied to the solution of the bound-state- 
like equations required by the LIT methods. However, 
due to differences in the asymptotics of the wave func- 
tions and in the strength distribution in the continuum 
achieved with the HO and HH expansions, the practi- 
cal implementation of the method, especially concerning 
the problems of convergence, might lead to difficulties. 
In particular, to ensure a small numerical uncertainty 
in the response function, obtained by numerical inver- 
sion |33 . l33 | , one has to achieve a very good accuracy in 
the calculation of the LIT. Consequently, it is necessary 
to find a range of HO frequencies, for which both the 
ground and the excited states of the system present good 
convergence properties. The actual choice of HQ depends 
on both the nucleus under consideration and the range of 
the transition operator. For "^He we find that frequencies 
in the range 12 MeV < hn < 28 MeV have the demanded 
characteristics for both the isovcctor dipole and isoscalar 
quadrupole excitations. 

Our main conclusion is that, in this benchmarking cal- 
culation with a semirealistic NN interaction, the NCSM 
is able to reach the level of precision required for the 
description of continuum responses via the LIT method 
and to obtain equivalent results as those found with the 
EIHH method. Because of the ability of NCSM to han- 
dle heavier mass nuclei than the EIHH, the present result 
opens the door for possible LIT investigations of heavier 
nuclei. However, the large model spaces needed in the 
NCSM calculations in order to achieve the necessary ac- 
curacy will require considerable thought and effort, but 
may be assisted by the use of effective field-theory two- 
and three-body interactions, which should be softer. 

Since LIT requires a good description of the ground- 
state wave function, and, at the same time, a fine 



enough discretization of the complex-energy contin- 
uum, one could use a mixed-mode approach in the 
spirit of Ref. (gJ. Thus, one could use basis states 
with a frequency that ensures fast convergence for 
the ground state, adding also basis states of smaller 
frequencies, which would allow a better discretization of 
the complex-energy continuum at low energies, and basis 
states corresponding to larger frequencies for a better 
description of the response at higher energies. This 
approach would require, however, a more involved nu- 
merical technique which can be adapted for our no-core 
calculations in a Jacobi basis. Moreover, such a method 
is not immediately applicable to Slater determinant basis 
codes, which are much more efficient for calculations in 
heavier nuclei, because of the more delicate treatment 
of the center-of-mass motion required in such a case. 
While such an approximation is under consideration, we 
leave its possible implementation for future work. 
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